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Abstract 

The recent detection of delayed, low energy emission from Gamma-Ray Burst (GRB) sources 
confirmed the cosmological origin of the bursts and provided support for models where GRBs are 
produced by the dissipation of the kinetic energy of relativistic fireballs. In this review, ultra- 
high-energy, > 10 19 eV, cosmic-ray and high energy, ~ 10 14 eV, neutrino production in GRBs 
is discussed in the light of recent GRB and cosmic-ray observations. Emphasis is put on model 
predictions that can be tested with operating and planned cosmic-ray and neutrino detectors. 
The predicted neutrino intensity, E^dN u /dE u ~ 3 x 10" 9 GeV/cm 2 s sr for 10 14 eV< E v < 
10 16 eV, implies that a km 2 neutrino detector would observe tens of events per year correlated 
with GRBs, and will be able to test for neutrino properties with an accuracy many orders 
of magnitude better than is currently possible. The predicted production rate of high-energy 
protons, which is consistent with that required to account for the observed ultra-high-energy 
cosmic-ray (UHECR) flux, implies that operating and planned cosmic-ray detectors can test the 
GRB model for UHECR production. If the predicted sources are found, cosmic-ray detectors 
will provide us with a technique to investigate the inter-galactic magnetic field. 
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1 Introduction 



The origin of Gamma- Ray Bursts (GRBs), bursts of 0.1 MeV — 1 MeV photons lasting for a few 
seconds, remained unknown for over 20 years, primarily because GRBs were not detected until 
the past year at wave-bands other than 7-rays (see [|[ for review of 7-ray observations). The 
isotropic distribution of bursts over the sky, revealed by observations of the BATSE detector on 
board the Compton Gamma-Ray Observatory, suggested that GRB sources lie at cosmological 
distances 0. Adopting a cosmological distance scale to GRB sources, general phenomenological 
considerations were used to argue that the bursts are produced by the dissipation of the kinetic 
energy of a relativistic expanding fireball (see [H] for reviews). 

The availability from the BeppoSAX satellite of accurate positions for GRBs shortly after 
their detection enabled during the past year the detection of GRB "afterglows," delayed X- 
ray [[|, optical |5j and radio || emission associated with GRB sources. Optical afterglow 
observations confirmed the cosmological origin of the bursts: Absorption lines detected in the 
afterglow of one burst set a lower limit z > 0.835 to its redshift and the redshifts of three 
GRB host-galaxies were determined, GRB970508 at z = 0.835 §, GRB980703 at z = 0.965 ||, 
and GRB971214 at z — 3.42 [[[(]]. The characteristics of observed GRB afterglows, the existence 
of which has been predicted by the fireball model |IT| , are broadly in agreement with fireball 
model predictions and therefore provide strong support for the model [12]. Furthermore, the 
detection || of predicted |13| afterglow radio scintillation directly demonstrates the relativistic 
expansion of the source of GRB970508 fll4"| . It should be noted, however, that despite the 
general success of the fireball model the underlying sources producing GRB fireballs remain 
unknown. 

Much like the underlying GRB source, the origin of ultra-high energy cosmic-rays (UHE- 
CRs), cosmic-rays of energy > 10 19 eV, is unknown (see for a recent review). Most of 
the sources of cosmic-rays that have been proposed have difficulties in accelerating particles 

which is in excess of 10 20 



17, 18, 19 



up to the highest observed energy [T6|, |i5] , which is in excess of 10 eV 
Furthermore, since the distance traveled by the highest energy particles must be smaller, due 
to interaction with radiation backgrounds p2], than 100 Mpc [231, their arrival directions are 



inconsistent with the position of any astrophysical object that is likely to produce high energy 
particles p4[ . Well before the recent confirmation of the hypothesis that GRBs are of cos- 
mological origin, it has been shown that cosmological fireballs are likely sources of > 10 19 eV 
protons. The physical conditions in the fireball dissipation region imply that protons may be 
Fermi accelerated in this region to energy > 10 20 eV |^5[ 26]. In addition, the average rate at 



which energy is emitted as 7-rays by GRBs is remarkably comparable to the energy generation 
rate of UHECRs in a model where UHECRs are produced by a cosmological distribution of 
sources [[25|, 27| . These two facts suggest that GRBs and UHECRs have a common originF]. The 
GRB fireball model for UHECR production makes several unique predictions |29], |30|, which can 
be tested with operating (HiRes [0TJ]), and planned (Auger [EJ, Telescope- Array UHECR 
detectors. Possibly the most interesting consequence of proton acceleration in GRB fireballs 
is the conversion of a significant fraction, ~ 10%, of the fireball energy to an accompanying 



2 Milgrom & Usov suggested |28 a GRB-UHECR association based on the overlap of the Fly's Eye highest 
energy cosmic ray arrival direction error box with the position error box of a bright GRB. A correlation between 
GRB and UHECR arrival directions is not expected, as explained in §4, for cosmological GRBs. 
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burst of ~ 10 eV neutrinos [j34|| - The predicted flux implies detection of tens of events per 
year correlated with GRBs in planned km 2 neutrino detectors (AMANDA-II and Deeplce |35|j 



extensions of the operating AMANDA detector 36 , ANTARES B7L NESTOR 



In this paper, high-energy cosmic-ray and neutrino production in GRBs is discussed in the 
light of recent GRB and ultra-high-energy cosmic-ray observations. In §2 the fireball model is 
briefly described, and implications of recent afterglow observations are discussed, which are of 
importance for high energy particle production. A more detailed discussion of GRBs and the 
fireball model is given in a separate contribution to these proceedings ||39|| . In section §3 recent 



UHECR observations, described in detail elsewhere in these proceedings p0| , are discussed 



The flux and spectrum measured by the Fly's Eye, Yakutsk, and AGASA experiments are 
compared with the prediction of a model where UHECRs are protons accelerated to high 
energy by Fermi shock acceleration in sources which are uniformly distributed in the universe. 
The main spectral feature predicted by such a model is a "GZK cutoff" P2"| , a suppression of 
UHECR flux above ~ 5 x 10 19 eV due to interaction of protons with the microwave background. 
We show that present data do not allow to confirm or rule out the existence of the predicted 
suppression. While both Fly's Eye and Yakutsk data show a suppression at 2a significance, a 
discrepancy may be emerging between this data and the results of the AGASA experiment: the 
Fly's Eye experiment and the Yakutsk experiment each report one event beyond 10 20 eV (close 
to the average number of 1.5 events predicted by the cosmological model), while the AGASA 
experiment reports 6 events for similar exposure. Much larger exposure than presently available 
is required to determine whether or not a GZK "cutoff" exists. 

The production in GRB fireballs of UHECRs is discussed in §4. Predictions of the GRB 
model for UHECR production, that can be tested with future UHECR experiments, are dis- 
cussed in §5. High energy neutrino production in fireballs and its implications for future high 
energy neutrino detectors are discussed in §6. The discussion in §4 §6 of UHECR and neutrino 
production in GRBs is similar to the analysis presented prior to the discovery of GRB afterglow 
[ p5| , |3"3|| . Some quantitative modifications are introduced due to the revised GRB energy scale, 
~ 10 53 erg (for isotropic emission) implied by afterglow observations compared to ~ 10 52 erg 
previously assumed, and due to some, yet inconclusive, evidence that the local GRB rate is 
lower than previously estimated, ~ 1 Gpc _3 yr _1 compared to ~ 10 Gpc~ 3 yr _1 . We also ad- 
dress some criticism of the GRB model of UHECR production recently made in the literature 



regarding the energy loss of protons escaping the fireball HO, and the acceleration process [42 



We show that this criticism is inapplicable to the model for UHECR production discussed here, 
which was proposed in [p5] . Finally, in §7 a summary is presented of the main points discussed 
in the paper. 



2 GRB fireballs and afterglow observations 
2.1 The fireball model 

General phenomeno logical considerations, based on 7-ray observations, indicate that, regardless 
of the nature of the underlying sources, GRBs are produced by the dissipation of the kinetic 
energy of a relativistic expanding fireball. The rapid rise time and short duration, ~ 1 ms, 
observed in some bursts HJ| imply that the sources are compact, with a linear scale r ~ 10 7 cm. 
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The high 7-ray luminosity required for cosmological bursts, L 1 ~ 10 52 erg s -1 , then results in 
an initially optically thick (to pair creation) plasma of photons, electrons and positrons, which 
expands and accelerates to relativistic velocities |14] . This is true whether the energy is released 



instantaneously, i.e. over a time scale ro/c, or as a wind over a duration comparable to the 
entire burst duration (~seconds). In fact, the hardness of the observed photon spectra, which 
extends to ~ 100 MeV, implies that the 7-ray emitting region must be moving relativistically, 
with a Lorentz factor V ~ 300 [fi"5|| , in order that the fireball pair-production optical depth be 



small for the observed high energy photons. 

If the observed radiation is due to photons escaping the fireball/wind as it becomes optically 
thin, two problems arise. First, the photon spectrum is quasi-thermal, in contrast with obser- 
vations. Second, the source size, ro ~ 10 7 cm, and the total energy emitted in gamma-rays, 
~ 10 53 erg, suggests that the underlying energy source is related to the gravitational collapse 
of ~ 1M Q object. Thus, the plasma is expected to be "loaded" with baryons which may be 
injected with the radiation or present in the atmosphere surrounding the source. A small bary- 
onic load, > 1O~ 8 M , increases the optical depth (due to Thomson scattering) so that most 
of the radiation energy is converted to kinetic energy of the relativistically expanding baryons 
before the plasma becomes optically thin [[|(| . To overcome both problems it was proposed |47 



that the observed burst is produced once the kinetic energy of the ultra-relativistic ejecta is 
re-randomized by some dissipation process at large radius, beyond the Thomson photosphere, 
and then radiated as 7-rays. Collision of the relativistic baryons with the inter-stellar medium 



47]], and internal collisions within the ejecta itself ||48|| , were proposed as possible dissipation 
processes. 

Most GRBs show variability on time scales much shorter than (typically one hundredth 
of) the total GRB duration. Such variability is hard to explain in models where the energy 
dissipation is due to external shocks |49]. [50[| . Thus, it is believed that internal collisions are 
responsible for the emission of gamma-rays. At small radius, the fireball bulk Lorentz factor, 
r, grows linearly with radius, until most of the wind energy is converted to kinetic energy and 
T saturates at V ~ 300. Variability of the source on time scale At, resulting in fluctuations 
in the wind bulk Lorentz factor V on similar time scale, then leads to internal shocks in the 
expanding fireball at a radius 

n « T 2 cAt = 3 x 10 13 r 2 00 At 10ms cm, (1) 

where T = 300r3oo, At = lOAtioms rns. If the Lorentz factor variability within the wind is 
significant, internal shocks would reconvert a substantial part of the kinetic energy to internal 
energy. It is assumed that this energy is then radiated as 7-rays by synchrotron and inverse- 
Compton emission of shock-accelerated electrons. For internal collisions, the observed gamma- 
ray variability time, ~ Vi/T 2 c ~ At, reflects the variability time of the underlying source, 
and the GRB duration reflects the duration over which energy is emitted from the source. 
Since the wind Lorentz factor is expected to fluctuate on time scales ranging from the shortest 
variability time At to the wind duration T, internal collisions will take place over a range of 
radii, r ~ = T 2 cAt to r ~ T 2 cT. A large fraction of bursts detected by BATSE show 
variability on the shortest resolved time scale, ~ 10 ms. Our choice in Eq. (0) of At = 10 ms as 
a representative value is therefore conservative, in the sense that the shortest variability time 
may be significantly smaller. In fact, recent analysis indicates that variability on ~ 1 ms is 
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common [|5lt . 



Internal shocks are expected to be "mildly" relativistic in the fireball rest frame ||48|| , i.e. 
characterized by Lorentz factor ji — 1 ~ 1 (since adjacent shells within the wind are expected 
to expand with Lorentz factors which do not differ by more than an order of magnitude). The 
internal shocks would therefore heat the protons to random velocities (in the wind frame) 7 P — 
1 ~ 1. The characteristic frequency of synchrotron emission is determined by the characteristic 
energy of the electrons and by the strength of the magnetic field. These are determined by 
assuming that the fraction of energy carried by electrons is £ e , implying a characteristic rest 
frame electron Lorentz factor j e = £ e (^pA«e), and that a fraction of the energy is carried by 
the magnetic field, implying 47rrfcT 2 B 2 /87T = where L is the total wind luminosity. Since 
the electron synchrotron cooling time is short compared to the wind expansion time, electrons 
lose their energy radiatively and L w L 7 /£ e . The characteristic observed energy of synchrotron 
photons, E 1 = Yh^leB / 'm e c, is therefore 



£ 7 ^0.1(W0.3) 1/2 (£ e /0.3) 



3/2 _ 



-1/2 
^7,52 



TsOO^lOms 



-MeV, 



(2) 



where L 7 = 10 52 L 7; 52erg/s. At present, there is no theory that allows the determination of the 
values of the equipartition fractions £ e and However, it is encouraging that for values close 
to equipartition, the photon energy predicted by the model is similar to that observed. 

As the fireball expands, it drives a relativistic shock (blastwave) into the surrounding gas, 
e.g. into the inter-stellar medium (ISM) gas if the explosion occurs within a galaxy. In what 
follows, we refer to the surrounding gas as "ISM gas," although the gas need not necessarily 
be inter-stellar. At early time, the fireball is little affected by the interaction with the ISM. 
At late time, most of the fireball energy is transferred to the ISM, and the flow approaches 
the self-similar blast-wave solution of Blandford & McKee H52] . At this stage a single shock 
propagates into the ISM, behind which the gas expands with Lorentz factor 
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(3) 



where E = 10 53 E 53 erg is the (isotropic) fireball energy, n = lni cm -3 is the ISM number 
density, and r = 10 17 r 17 cm is the shell radius. The expansion becomes self-similar at a radius 
r where two conditions are met: the Lorentz factor Tbm^) inferred from the self-similar solution 
is smaller than the initial Lorentz factor T, and the width of the shell into which the shocked 
ISM is compressed in the self-similar solution, « r /10T 2 BM (r), is larger than the initial fireball 
shell width cT. The first and second conditions are met for r > rp = (17 E /\QTxY 2 nm p c 2 ) 1 ^ 



and r > tt = (10 x 17EcT/167mm p c 
external-shock flow occurs at 



2U/4 



respectively. Thus, the transition to self-similar, 



r e = 5.2 x 10 
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E< 



i/i 



max 



l.OTi 74 , 1.2 



Ei 



1/12 



,-2/3 
300 



cm, 



(4) 



where T = 17\ s. 

Internal, mildly-relativistic shocks within the fireball shell result both from variability of 
the source, at < r < r e , and from the interaction of the fireball with the surrounding gas, 
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during the transition to a self-similar expansion at r ~ r e , where reverse shocks propagate 
into the expanding fireball ejecta and decelerate it. From Eq. (|j) we infer that for typical 
fireball parameters the two conditions required for transition to self-similar external shock flow 
are satisfied at a similar radius, r e « r-p « tt- This implies that significant deceleration of 
the fireball shell does not take place prior to the transition to self-similar behavior. This, in 
turn, implies that the reverse shocks propagating into fireball ejecta are only mildly relativistic. 
Thus, the characteristics of internal shocks due to interaction with surrounding gas are similar 
to those of internal shocks due to variability on time scale ~ T. In the discussion that follows 
we therefore do not discuss the reverse shocks separately from the internal shocks. 

The shock driven into the ISM continuously heats new gas, and produces relativistic elec- 
trons that may produce the delayed radiation, "afterglow" , observed on time scales of days to 
months. As the shock-wave decelerates, the emission shifts to lower frequency with time. Since 
proton acceleration to high energy takes place only in the internal shocks, r < r e (see §4), we 
do not discuss further the theory of afterglow emission. 



2.2 Afterglow observations 

Afterglow observations confirmed, as discussed in the Introduction, the cosmological origin of 
GRBs, and are consistent with delayed GRB emission being synchrotron radiation of elec- 
trons accelerated to high energy in the highly relativistic shock driven by the fireball into its 
surrounding gas. Since we are interested mainly in the earlier, internal collision phase, of 
fireball evolution, we do not discuss afterglow observations in detail. We note, however, two 
implications of afterglow observations which are of importance for the discussion of UHECR 
production. 

The first implication is related to the GRB energy scale. The gamma-ray energy emitted 
by the three GRBs with measured redshifts in the energy range of 20 keV to 2 MeV (assuming 
spherical symmetry) is ss 10 52 erg, ps 0.8 x 10 53 erg and w 3 x 10 53 erg for GRB970508 
(z = 0.835), GRB980703 (z = 0.966) and GRB971214 (z = 3.42) respectively (here, and 
throughout the paper, we assume an open universe, Q = 0.2, A = 0, and H = 75 km/s Mpc). 
This implies that GRBs are not "standard candles," and that the characteristic gamma-ray 
energy (luminosity) is ~ 10 53 erg (~ 10 52 erg/s) rather than ~ 10 52 erg (~ 10 51 erg/s) as 
commonly assumed in the past. Performing a detailed analysis Mao & Mo (1998) find, for 
example, that the typical luminosity of observed GRBs is ~ 10 52 erg/s (Note that Mao & Mo 
use Hq = 100 km/Mpc s and quote luminosities and energies in the 50-300 keV band only). 

The second implication relates to the GRB rate. Krumholtz, Thorsett & Harrison (see 



also [5"4"|) have demonstrated that, based on the data, it is impossible to distinguish between 
models where the GRB rate per unit comoving volume is independent of redshift, and models 
where it evolves rapidly, e.g. following star formation rate (previous claims to the contrary were 
based on the assumption, now known to be invalid, that GRBs are standard candles). Most 



observed GRBs originate at the redshift range z ~ 1 to z ~ 2 [[5^, [54[], and the observed GRB 
rate essentially determines the GRB rate per unit volume at that redshift. The present rate is 
less well constrained and ranges from -Rgrb ~ 1/Gpc 3 yr, assuming the GRB rate evolves rapidly 
as the star- formation rate, to -Rgrb ~ 10/Gpc 3 yr, assuming the GRB rate is independent of 
redshift. There is some evidence supporting the hypothesis that the GRB rate follows the star 
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Figure 1: The UHECR flux expected in a cosmological model, where high-energy protons are 
produced at a rate E 2 dficR/dE = 0.8 x 10 44 erg/Mpc 3 yr, compared to the Fly's Eye, Yakutsk 
and AGASA data, la flux error bars are shown. The highest energy points are derived assuming 
the detected events (1 for Fly's Eye and Yakutsk, 4 for AGASA) represent a uniform flux over 
the energy range 10 20 eV-3 x 10 20 eV. 



formation rate |55fl . The evidence is, however, not yet conclusive. 



3 UHECR observations and their implications 



Fly's Eye |J9] and AGASA p0, |21| results confirm the flattening of the cosmic-ray spectrum 
at ~ 10 19 eV, evidence for which existed in previous experiments with weaker statistics fI7 |. 
Fly's Eye data is well fitted in the energy range 10 17,6 eV to 10 19 ' 6 eV by a sum of two power 
laws: A steeper component, with differential number spectrum J oc E~ 3,50 , dominating at lower 
energy, and a shallower component, J oc E~ 2S1 , dominating at higher energy, E > 10 19 eV. 
The flattening of the spectrum, combined with the lack of anisotropy and the evidence for a 
change in composition from heavy nuclei at low energy to light nuclei (protons) at high energy 



H], |TJ|] , suggest that an extra-Galactic source of protons dominates the flux at high energy. 

In Fig. 1 we compare the UHECR spectrum, reported by the Fly's Eye, the Yakutsk, and 
the AGASA experiments |T5| , [T8] , with that expected from a homogeneous cosmological 
distribution of sources, each generating a power law differential spectrum of high energy protons 
dN/dE oc E~ 2 . This spectrum is expected for Fermi shock acceleration (see §4.1 for 
discussion of the GRB UHECR model). The absolute flux measured at 3 x 10 18 eV differs 
between the various experiments, corresponding to a systematic ~ 10% (~ 20%) over-estimate 
of event energies in the AGASA (Yakutsk) experiment compared to the Fly's Eye experiment 
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(see also [20]). In Fig. 1, the Yakutsk energy normalization is used. For the model calculation, 
an open universe, Q = 0.2, A = and H = 75km/Mpc s were assumed. The calculation 
is similar to that described in |p7| |. The generation rate of cosmic-rays (per unit comoving 
volume) was assumed to evolve rapidly with redshift following the luminosity density evolution 
of QSOs |57| ) which is also similar to that describing the evolution of star formation rate [58| : 



ticr(z) oc (1 + z) a with a ~ 3 [^fj] at low redshift, z < 1.9, ncn(z) = Const, for 1.9 < z < 2.7, 
and an exponential decay at z > 2.7 [60|. The cosmic-ray spectrum at energy > 10 19 eV is 



little affected by modifications of the cosmological parameters or of the redshift evolution of 
cosmic-ray generation rate. This is due to the fact that cosmic-rays at this energy originate 
from distances shorter than several hundred Mpc. The spectrum and flux at E > 10 19 eV is 
mainly determined by the present (z = 0) generation rate and spectrum, which in the model 
shown in Fig. 1 is E 2 (dhcR/dE) z=0 = 0.8 x 10 44 erg/Mpc 3 yr. 

The suppression of model flux above ~ 10 19 ' 7 eV, compared to a power-law extrapolation 
of the flux at low energy, is due to energy loss of high energy protons in interaction with the 
microwave background, i.e. to the "GZK cutoff" |[22|| . This is the characteristic signature 
of cosmological source distribution. It is clear from Fig. 1 that present data does not allow 
to confirm or rule out the existence of the "cutoff" with high confidence. Nevertheless, some 
evidence for the cutoff does exit. Both Fly's Eye and Yakutsk data show a deficit in the number 
of events detected above 10 19 7 eV, compared to the number expected based on extrapolation 
of the Fly's Eye shallower power- law fit for E < 10 19 ' 6 eV, J oc E~ 2 - 61 . The deficit is, however, 
only at a 2cr confidence level |[25|| . The AGASA data is consistent with Fly's Eye and Yakutsk 
results below 10 20 eV. A discrepancy may be emerging, see Fig. 1, at higher energy, > 10 20 eV, 
where the Fly's Eye and Yakutsk experiments detect 1 event each, and the AGASA experiment 
detects 6 events for similar exposure (a ~ 2a discrepancy). 

We therefore conclude, that a scenario where UHECRs are produced by a cosmological distri- 
bution of sources, generating high energy protons at a rate E 2 (dh C R/dE) z=0 m 10 44 erg/Mpc 3 yr, 
is consistent with the observed flux and spectrum of cosmic-rays in the energy range of 10 19 eV 
to 10 20 eV. The flux predicted by this model above 10 20 eV is consistent with that measured by 
the Fly's Eye and the Yakutsk experiments, while AGASA results suggest a higher flux at this 
energy. The statistical significance of the discrepancy between the experiments (or between 
the AGASA results and model prediction above 10 20 eV) is not high. Clearly, much larger 
exposure than presently available is required to accurately determine the UHECR spectrum 
and flux above 5 x 10 19 eV. 



4 UHECRs from GRB fireballs 

4.1 Fermi acceleration in GRBs 

In the fireball model, the observed radiation is produced, both during the GRB and the af- 
terglow, by synchrotron emission of shock accelerated electrons. In the region where electrons 
are accelerated, protons are also expected to be shock accelerated. This is similar to what is 
thought to occur in supernovae remnant shocks, where synchrotron radiation of accelerated 
electrons is the likely source of non-thermal X-rays (recent ASCA observations give evidence 
for acceleration of electrons in the remnant of SN1006 to 10 14 eV ^TJ), and where shock accel- 
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eration of protons is believed to produce cosmic rays with energy extending to ~ 10 15 eV 



e.g., |62| for review). Thus, it is likely that protons, as well as electrons, are accelerated to high 
energy within GRB fireballs. Let us consider the constraints that should be satisfied by the 
fireball parameters in order to allow acceleration of protons to ~ 10 20 eV. 

We consider proton Fermi acceleration in fireball internal shocks, which take place as the 
fireball expands over a range of radii r« ~ 10 14 cm < r < r e ~ 10 16 cm [cf. Eqs. (2,1])]. As 
mentioned in §2.1, internal shocks are due to variability of the source, on time scales ranging 
from the shortest variability time At ~ 10 ms to the wind duration T ~ 10 s. In addition, 
internal shocks arise also due to interaction with ambient gas. The characteristics of the latter 
internal shocks, which occur at r ~ r e during the transition to self-similar expansion, are similar 
to those of internal collisions due to variability on time scale T [see discussion following Eq. (fj)]. 
Internal shocks are, in the wind rest-frame, "mildly relativistic," i.e. characterized by Lorentz 
factors 7« — 1 ~ 1. We therefore expect results related to particle acceleration in sub-relativistic 
shocks to be valid for the present scenario. The most restrictive requirement, which rules out 
the possibility of accelerating particles to energy ~ 10 20 eV in most astrophysical objects, is 
that the particle Larmor radius Rl should be smaller than the system size In our scenario 
we must apply a more stringent requirement, namely that Rl should be smaller than the largest 
scale / over which the magnetic field fluctuates, since otherwise Fermi acceleration may not be 
efficient. We may estimate / as follows. The comoving time, i.e. the time measured in the 
wind rest frame, is t = r/Tc. Thus, regions separated by a comoving distance larger than r/T 
are causally disconnected, and the wind properties fluctuate over comoving length scales up to 
/ ~ r/T. We must therefore require Rl < r/T. A somewhat more stringent requirement is 
related to the wind expansion. Due to expansion the internal energy is decreasing and therefore 
available for proton acceleration (as well as for 7-ray production) only over a comoving time 
t ~ r/Tc. The typical Fermi acceleration time is t a = f'RL/c/3 2 , where f3c is the Alfven velocity 
and / ~ 1 |63|, In our scenario (3 ~ 1 leading to the requirement JRl < r/T. This condition 
sets a lower limit to the required comoving magnetic field strength. This limit may be stated 
as a radius independent lower limit to the ratio of magnetic field and electron energy densities 



25 



W£e > 0.02/ 2 r 2 00 £ 2 oL-i, (5) 



where E p = 10 20 -E Pi 2o eV is the accelerated proton energy. 

The accelerated proton energy is also limited by energy loss due to synchrotron radiation and 
interaction with fireball photons. As discussed in §6, the dominant energy loss process is syn- 
chrotron cooling. The condition that the synchrotron loss time, t sy = (67rmpC 3 /aTfnl)E^ 1 B^ 2 , 
should be smaller than the acceleration time sets an upper limit to the magnetic field strength. 
Since the equipartition field decreases with radius, B ep . oc r~ 2 , the upper limit on the magnetic 
field may be satisfied simultaneously with (|5|) provided that the internal collisions occur at large 



enough radius [25 



r > 10 12 / 2 r 3rj 2 E p 3 20 cm. (6) 
Since collisions occur at radius r m T 2 cAt, the condition (|6p is equivalent to a lower limit on F 

T > 130/V^Airol 4 . (7) 

^From Eqs. (0) and (|7p, we infer that a dissipative ultra-relativistic wind, with luminosity 
and variability time implied by GRB observations, satisfies the constraints necessary to allow 
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the acceleration of protons to energy > 10 20 eV, provided that the wind bulk Lorentz factor is 
large enough, F > 100, and that the magnetic field is close to equipartition with electrons. The 
former condition, T > 100, is remarkably similar to that inferred based on the 7-ray spectrum, 
and r ~ 300 is the "canonical" value assumed in the fireball model. The latter condition, 
magnetic field close to equipartition, is commonly assumed to be valid in order to account for 
the observed 7-ray emission [see, e.g., Eq. @]. 

Finally, two points should be clarified. First, it has recently been claimed that ultra-high 
energy protons would lose most of their energy adiabatically, i.e. due to expansion, before 
they escape the fireball This claim is based on the assumptions that internal shocks, and 
therefore proton acceleration, occur at r ~ only, and that subsequently, r« < r < r e , the 
fireball expands adiabatically. Under these assumptions, protons would lose most their energy 
by the time they escape, at r ~ r e . However, as emphasized both in this section and in §2.1, 
internal shocks are expected to occur at all radii r, < r < r e , and in particular at r ~ r e 
during the transition to self-similar expansion. Thus, proton acceleration to ultra-high energy 



is expected to operate at all radii up to r 



where ultra-high energy particles escape. 



Second, it has recently been pointed out in ||2| that the conditions at the external shock 
driven by the fireball into the ambient gas are not likely to allow proton acceleration to ultra- 
high energy. Although correct, this observation is irrelevant for the acceleration in internal 



shocks, the scenario considered for UHECR production in GRBs in both [25] and [26 



4.2 UHECR flux and spectrum 

We have shown in §3, that the present rate at which energy should be produced as > 10 19 eV pro- 
tons by cosmological cosmic-ray sources in order to produce the observed flux is E 2 {dhcnl dE) z=Q 
10 44 erg Mpc _3 yr _1 . This rate is, remarkably, comparable to that produced in 7-rays by cos- 
mological GRBs. The typical GRB 7-ray energy, E ~ 10 53 erg, and the present [z = 0) GRB 
rate, which is estimated to be in the range of Rgrb = 1/Gpc 3 yr to Rgrb = 10/Gpc 3 yr (see 
discussion in §2.2), implies a present energy generation rate in the range ~ 10 44 erg Mpc" 3 yr" 1 
to ~ 10 45 erg Mpc _3 yr _1 . In addition, since protons are accelerated in the GRB model to high 
energy by internal shocks, which in the fireball frame are sub-relativistic, we may expect a 
generation spectrum dN/dE oc E~ 2 , consistent with UHECR observations (see §3.1). 

Thus, GRB fireballs would produce UHECR flux and spectrum consistent with that ob- 
served, provided the efficiency with which the wind kinetic energy is converted to 7-rays, and 
therefore to electron energy, is similar to the efficiency with which it is converted to proton en- 
ergy, i.e. to UHECRs [2S|. There is, however, one additional point which requires consideration 



£2). The energy of the most energetic cosmic ray detected by the Fly's Eye experiment is in 
excess of 2 x 10 20 eV, and that of the most energetic AGASA event is ~ 2 x 10 20 eV. On a cos- 
mological scale, the distance traveled by such energetic particles is small: < lOOMpc (50Mpc) 
for the AGASA (Fly's Eye) event (e.g., f23|). Thus, the detection of these events over a ~ 5yr 



period can be reconciled with the rate of nearby GRBs, ~ 1 per 100 yr out to lOOMpc, only if 
there is a large dispersion, > lOOyr, in the arrival time of protons produced in a single burst 
(This implies that if a direct correlation between high energy CR events and GRBs, as recently 
suggested in p8| , is observed on a ~ lOyr time scale, it would be strong evidence against a 



cosmological GRB origin of UHECRs). 
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The required dispersion is likely to occur due to the combined effects of deflection by random 
magnetic fields and energy dispersion of the particles |[25|| . Consider a proton of energy E 
propagating through a magnetic field of strength B and correlation length A. As it travels 
a distance A, the proton is typically deflected by an angle a ~ \/Rl, where Rl = E/eB 
is the Larmor radius. The typical deflection angle for propagation over a distance D is 9 S ~ 
{D/\) l/2 \/R L . This deflection results in a time delay, compared to propagation along a straight 
line, 

r{E, D) « 9 2 D/Ac « 2 x 10 5 £ 2 - 2 I^ 00 Ampc£*g 5*, (8) 

where D = 100-DiooMpc, A = 1Am pc Mpc and B = 10~ 9 £> n G G. Here, we have chosen numerical 
values corresponding to the current upper bound on the inter-galactic magnetic field, BX 1 ^ 2 < 
]. The random energy loss UHECRs suffer as they propagate, owing to the 
production of pions, implies that at any distance from the observer there is some finite spread in 
the energies of UHECRs that are observed with a given fixed energy. For protons with energies 
> 10 20 eV the fractional RMS energy spread is of order unity over propagation distances in the 
range 10 — lOOMpc (e.g. [p3fl ). Since the time delay is sensitive to the particle energy, this 
implies that the spread in arrival time of UHECRs with given observed energy is comparable 
to the average time delay at that energy t(E, D) (This result has been confirmed by numerical 
calculations in p5|j ). Thus, the required time spread, r > 100 yr, is consistent with the upper 
bound, r < 2 x 10 5 yr, implied by the present upper bound to the inter-galactic magnetic field. 



io- 9 g Mpc 1 ^ m 



5 GRB model predictions for UHECR experiments 
5.1 The Number and Spectra of Bright Sources 

The initial proton energy, necessary to have an observed energy E, increases with source dis- 
tance due to propagation energy losses. The rapid increase of the initial energy after it exceeds, 
due to electron-positron production, the threshold for pion production effectively introduces 
a cutoff distance, D C (E), beyond which sources do not contribute to the flux above E. The 



function D C (E) is shown in Fig. 2 (taken from [29]). Since D C (E) is a decreasing function of 
E, for a given number density of sources there is a critical energy E c , above which only one 
source (on average) contributes to the flux. In the GRB model E c depends on the product of 
the burst rate Rgrb and the time delay. The number of sources contributing, on average, to 
the flux at energy E is [f29f| 

47T 

N(E) = —R GRB D c (E) 3 r [E, D C (E)} , (9) 
5 

and the average intensity resulting from all sources is 

J(E) = ^R GRB ^D C (E) , (10) 

where dn p /dE is the number per unit energy of protons produced on average by a single burst 
(this is the formal definition of D C (E)). The critical energy E c is given by 

47T 

—R grb D c (E c ) 3 t[E c ,D c (E c )] = 1 . (11) 
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Figure 2: Results of a Monte-Carlo realization of the bursting sources model, with E c = 
1.4 x 10 20 eV: Thick solid line- overall spectrum in the realization; Thin solid line- average 
spectrum, this curve also gives D C (E); Dotted lines- spectra of brightest sources at different 
energies. 



E c , the energy beyond which a single source contributes on average to the flux, depends 
on the unknown properties of the intergalactic magnetic field, r oc B 2 X. However, the rapid 
decrease of D C (E) with energy near 10 20 eV, see Fig. 2, implies that E c is only weakly dependent 
on the value of B 2 X. In The GRB model, the product R grb t(D = lOOMpc, E = 10 20 eV) is 
approximately limited to the range 10~ 6 Mpc~ 3 to 10~ 3 Mpc -3 (The lower limit is set by the 
requirement that at least a few GRB sources be present at D < 100 Mpc, and the upper limit 
by the Faraday rotation bound 5A 1/2 < 10" 9 G Mpc 1/2 
corresponding range of values of E c is 10 20 eV < E c < 3 x 10 20 eV. 

Fig. 2 presents the flux obtained in one realization of a Monte-Carlo simulation described 
by Miralda-Escude & Waxman of the total number of UHECRs received from GRBs at 
some fixed time. For each realization the positions (distances from Earth) and times at which 
cosmological GRBs occurred were randomly drawn, assuming, an intrinsic proton generation 
spectrum dN p /dE p oc E~ 2 , and E c = 1.4 x 10 20 eV. Most of the realizations gave an overall 
spectrum similar to that obtained in the realization of Fig. ^ when the brightest source of 
this realization (dominating at 10 20 eV) is not included. At E < E c , the number of sources 
contributing to the flux is very large, and the overall UHECR flux received at any given time 
is near the average (the average flux is that obtained when the UHECR emissivity is spatially 
uniform and time independent). At E > E c , the flux will generally be much lower than the 
average, because there will be no burst within a distance D C (E) having taken place sufficiently 
recently. There is, however, a significant probability to observe one source with a flux higher 



and R GRB < 10/ GpC». The 
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than the average. A source similar to the brightest one in Fig. |2| appears ~ 5% of the time. 

At any fixed time a given burst is observed in UHECRs only over a narrow range of energy, 
because if a burst is currently observed at some energy E then UHECRs of much lower energy 
from this burst have not yet arrived, while higher energy UHECRs reached us mostly in the 
past. As mentioned above, for energies above the pion production threshold, £~5x 10 19 eV, 
the dispersion in arrival times of UHECRs with fixed observed energy is comparable to the 
average delay at that energy. This implies that the spectral width AE of the source at a given 
time is of order the average observed energy, AE ~ E. Thus, bursting UHECR sources should 
have narrowly peaked energy spectra, and the brightest sources should be different at different 
energies. For steady state sources, on the other hand, the brightest source at high energies 
should also be the brightest one at low energies, its fractional contribution to the overall flux 
decreasing to low energy only as Dc(E)^ 1 . A detailed numerical analysis of the time dependent 
energy spectrum of bursting sources is given in K| B7 . 



5.2 Spectra of Sources at E < 4 x 10 19 eV 

The detection of UHECRs above 10 20 eV imply that the brightest sources must lie at distances 
smaller than lOOMpc. UHECRs with E < 4 x 10 19 eV from such bright sources will suffer energy 
loss only by pair production, because at E < 5 x 10 19 eV the mean-free-path for pion production 
interaction (in which the fractional energy loss is ~ 10%) is larger than lGpc. Furthermore, 
the energy loss due to pair production over lOOMpc propagation is only ~ 5%. 

In the case where the typical displacement of the UHECRs due to deflections by inter- 
galactic magnetic fields is much smaller than the correlation length, A ^> D9 S (D, E) ~ D(D / A) 1 / 2 
all the UHECRs that arrive at the observer are essentially deflected by the same magnetic field 
structures, and the absence of random energy loss during propagation implies that all rays with 
a fixed observed energy would reach the observer with exactly the same direction and time 
delay. At a fixed time, therefore, the source would appear mono-energetic and point-like. In 
reality, energy loss due to pair production results in a finite but small spectral and angular 
width, AE/E ~ 56 /6 a < 1% 0. 

In the case where the typical displacement of the UHECRs is much larger than the corre- 
lation length, A <C D6 S (D, E), the deflection of different UHECRs arriving at the observer are 
essentially independent. Even in the absence of any energy loss there are many paths from the 
source to the observer for UHECRs of fixed energy E that are emitted from the source at an 
angle 9 < 9 S relative to the source-observer line of sight. Along each of the paths, UHECRs 
are deflected by independent magnetic field structures. Thus, the source angular size would be 
of order 9 S and the spread in arrival times would be comparable to the characteristic delay r, 
leading to AE/E ~ 1 even when there are no random energy losses. The observed spectral 
shape of a nearby (D < lOOMpc) bursting source of UHECRs at E < 4 x 10 19 eV was derived 
for the case A <C D9 S (D, E) in |]3*0fl , and is given by 



dN °° 

^cx£(-l) n+1 ™ 2 exp 



2n 2 n 2 E 2 



(12) 



where E (t,D) = De(2B 2 \/3ct)^ 2 . For this spectrum, the ratio of the RMS UHECR energy 
spread to the average energy is 30% 
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B[G] 



Figure 3: The line 9 S D = A for a source at 30Mpc distance observed at energy E ~ 10 19 eV 
(dot-dash line), shown with the Faraday rotation upper limit BX 1 / 2 < 1(T 9 G Mpc 1/2 (solid 
line), and with the lower limit BX 1 ^ 2 > 10~ n G Mpc 1,/2 required in the GRB model. 

Fig. |3| shows the line 9 S D = X in the B — A plane, for a source at a distance D = 30Mpc 
observed at energy E ~ 10 19 eV. Since the 6 S D = X line divides the allowed region in the plane 
at A ~ IMpc, measuring the spectral width of bright sources would allow to determine if the 
field correlation length is much larger, much smaller, or comparable to IMpc. 



6 High energy Neutrinos 

6.1 Neutrino production 

6.1.1 Neutrinos at energies ~ 10 14 eV 

Protons accelerated in the fireball to high energy lose energy through photo-meson interaction 
with fireball photons. The decay of charged pions produced in this interaction, tt + —>■ [i + + v ll — ► 
e + + v e + 17^ + z/ M , results in the production of high energy neutrinos. The neutrino spectrum is 
determined by the observed gamma-ray spectrum, which is well described by a broken power- 
law, dN^/dE y oc E~@ with different values of (3 at low and high energy |68|| . The observed break 
energy (where (3 changes) is typically E h ~ lMeV, with (3 ~ 1 at energies below the break and 
(3 ~ 2 above the break. The interaction of protons accelerated to a power-law distribution, 
dN p /dE p oc E~ 2 , with GRB photons results in a broken power law neutrino spectrum |34 



dN u /dE u oc E~P with (3 = 1 for E v < E%, and (3 = 2 for E v > E h v . The neutrino break energy 
E h v is fixed by the threshold energy of protons for photo-production in interaction with the 
dominant ~ 1 MeV photons in the GRB, 

£>5x 10 14 r 2 00 (^/lMeV)- 1 eV. (13) 

The normalization of the flux is determined by the efficiency of pion production. As shown 



in 



p4[ , the fraction of energy lost to pion production by protons producing the neutrino flux 



14 



above the break, is essentially independent of energy and is given by 

f ~ n 2 Ll,b2 (14) 

Jn • (^/lMeV)r| 00 At 10ms - 1 ) 

Thus, acceleration of protons to high energy in internal fireball shocks would lead to conversion 
of a significant fraction of proton energy to high energy neutrinos. 

If GRBs are the sources of UHECRS, then using Eq. ([TJ]) and the UHECR generation rate 
implied by observations, the expected GRB neutrino flux is 

E 2 § m E 2 &- ps E 2 § 

~ 1.5 x 1(T 9 (J^J minl^^/^jGeVcm^V^r" 1 . (15) 



The GRB neutrino flux can be estimated directly from the observed gamma-ray fluence. 
The Burst and Transient Source Experiment (BATSE) measures the GRB fluence _F 7 over 
two decades of photon energy, ~ 0.02MeV to ~ 2MeV, corresponding to a decade of radiating 
electron energy (the electron synchrotron frequency is proportional to the square of the electron 
Lorentz factor). If electrons carry a fraction f e of the energy carried by protons, then the muon 
neutrino fluence of a single burst is E 2 dN v j dE v ps 0.25(/ 7r // e )F 7 / ln(10). The average neutrino 
flux per unit time and solid angle is obtained by multiplying the single burst fluence with the 
GRB rate per solid angle, ~ 10 3 bursts per year over An sr. Using the average burst fluence 
F 7 = 10~ 5 erg/cm 2 , we obtain a muon neutrino flux E 2 $> u ps 3 x 10 _9 (/^// e )GeV/cm 2 ssr. Thus, 
the neutrino flux estimated directly from the gamma-ray fluence agrees with the estimate fllBD 
based on the cosmic-ray production rate. 

6.1.2 Neutrinos at high energy > 10 16 eV 

The neutrino spectrum ( fllf) is modified at high energy, where neutrinos are produced by the 
decay of muons and pions whose life time r Mj7r exceeds the characteristic time for energy loss 



due to adiabatic expansion and synchrotron emission |34|, [70| |69|. The synchrotron loss time 
is determined by the energy density of the magnetic field in the wind rest frame. For the 
characteristic parameters of a GRB wind, the muon energy for which the adiabatic energy loss 
time equals the muon life time, E^, is comparable to the energy E s at which the life time equals 
the synchrotron loss time, r*. For pions, E° > E^. This, and the fact that the adiabatic loss 
time is independent of energy and the synchrotron loss time is inversely proportional to energy, 
imply that synchrotron losses are the dominant effect suppressing the flux at high energy. The 
energy above which synchrotron losses suppress the neutrino flux is 

E ^,"e) ^ f t r /M-1/2P2 A, , vb H A A „ f 10 > for *V V e 



El 



(e^ 7 ,52/e e )- 1/2 r 2 oAtio ms (^/ 1MeV ) >< { ioo, for (16) 



We note, that the results presented above were derived using the "A-approximation," i.e. 
assuming that photo-meson interactions are dominated by the contribution of the A-resonance. 
It has recently been shown [|7T|, that for photon spectra harder than dN^/dE^ oc E~ 2 , the 
contribution of non-resonant interactions may be important. Since in order to interact with the 
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hard part of the photon spectrum, E 1 < E^, the proton energy must exceed the energy at which 
neutrinos of energy E h v are produced, significant modification of the A- approximation results 
is expected only for E u ^> E b y) where the neutrino flux is strongly suppressed by synchrotron 
losses. 

So far, we have discussed neutrino production in internal shocks due to variability on the 
shortest time scale, At ~ 10 ms. Internal collisions due to variability on longer time scales, 
At < 5t < T, are less efficient in producing neutrinos, f n oc St^ 1 [cf. Eq. (pBD1, since the 
radiation energy density is lower at larger collision radii. However, at larger radii synchrotron 
losses cut off the spectrum at higher energy, E s (5t) oc St [cf. Eq. ( fL6|) l. Collisions at large radii 
therefore result in extension of the neutrino spectrum of Eq. ( |15|) to higher energy, beyond the 
cutoff energy Eq. ([16]), 

El$ v ^E-\ E U >E S U . (17) 

The neutrino flux from GRBs is small above 10 19 eV, and a neutrino flux comparable to 
the 7-ray flux is expected only below ~ 10 17 eV, in agreement with the results of ref. 



Our result is not in agreement, however, with that of ref. [72], where a much higher flux at 



10 19 eV is obtained based on the equations of ref. H34j| , which are the same equations as used 



hereQ. Finally, we note that, contrary to the claim in |70], there is no contradiction between 
production of high-energy protons above ~ 3 x 10 20 eV and a break in the neutrino spectrum 
at ~ 10 16 eV [cf. Eqs. (|g[T§)]. 

6.2 Implications 

The high energy neutrinos predicted in the dissipative wind model of GRBs may be observed by 
detecting the Cerenkov light emitted by high energy muons produced by neutrino interactions 



below a detector on the surface of the Earth (see [73[ for a recent review). The probability P Ufl 
that a neutrino would produce a high energy muon in the detector is approximately given by 
the ratio of the high energy muon range to the neutrino mean free path. At the high energy we 
are considering, P Vfl ~ 10 _6 (ej,/lTeV) |7!|. Using Eq. ([T5|), this implies a detection rate of ~ 20 
neutrino induced muons per year for a 1km 2 detector (over An sr). As discussed in |34 |, one 



may look for neutrino events in angular coincidence, on degree scale, and temporal coincidence, 



on time scale of seconds, with GRBs. Several authors |74[ have recently emphasized the effect 
on muon detection rate of neutrino absorption in the Earth. This effect is not large for GRB 
neutrinos, since most of the signal comes from neutrinos of energy ~ 10 14 eV. At this energy, 
the flux of upward moving muons is reduced due to absorption by 36% J7i[ , and the total 
(in sr) flux is reduced by only 18%. 

Detection of neutrinos from GRBs could be used to test the simultaneity of neutrino and 
photon arrival to an accuracy of ~ 1 s (~ 1 ms for short bursts), checking the assumption 
of special relativity that photons and neutrinos have the same limiting speed [The time delay 
for neutrino of energy 10 14 eV with mass m v traveling 100 Mpc is only ~ 10 _11 (m I/ /10 eV) 2 s]. 
These observations would also test the weak equivalence principle, according to which photons 

3 The parameters chosen in |72| are L 7 ~ 10 50 erg/s, At = 10s, and T = 100. Using equation (4) of ref. fl34f , 
which is the same as Eq. ( |l4|) of the present paper, we obtain for these parameters /, = 1.6 x 10~ 4 , while the 
author of [|72| obtains, using the same equation, = 0.03. 
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and neutrinos should suffer the same time delay as they pass through a gravitational potential. 
With 1 s accuracy, a burst at 100 Mpc would reveal a fractional difference in limiting speed 
of 10~ 16 , and a fractional difference in gravitational time delay of order 10~ 6 (considering the 
Galactic potential alone). Previous applications of these ideas to supernova 1987A (see for 
review), where simultaneity could be checked only to an accuracy of order several hours, yielded 
much weaker upper limits: of order 10~ 8 and 10~ 2 for fractional differences in the limiting speed 
f?6|1 and time delay |77| respectively. 

The model discussed above predicts the production of high energy muon and electron neu- 
trinos with a 2:1 ratio. If vacuum neutrino oscillations occur in nature, then neutrinos that get 
here should be almost equally distributed between flavors for which the mixing is strong. In 
fact, if the atmospheric neutrino anomaly has the explanation it is usually given, oscillation to 



z/ T 's with mass ~ 0.1 eV Jf8|, then one should detect equal numbers of z/ M 's and u^s. Upgoing 
r's, rather than /i's, would be a distinctive signature of such oscillations. Since i/ r 's are not 
expected to be produced in the fireball, looking for r's would be an "appearance experiment" 
(z/ r 's may be produced by photo-production of charmed mesons; However, the high photon 
threshold, ~ 50GeV, and low cross-section, ~ 1/ib JF{|, for such reactions imply that the ratio 



of charmed meson to pion production is ~ 10 ). To allow flavor change, the difference in 
squared neutrino masses, Am 2 , should exceed a minimum value proportional to the ratio of 
source distance and neutrino energy ]75]. A burst at 100 Mpc producing 10 14 eV neutrinos can 
test for Am 2 > 10~ 16 eV 2 , 5 orders of magnitude more sensitive than solar neutrinos. Note, 
that due to the finite pion life time, flavor mixing would be caused by de-coherence, rather than 
by real oscillations, for neutrinos with masses > O.leV. 



7 Summary 

Afterglow observations confirmed the cosmological origin of GRBs and provide support for the 
fireball model (§1 ,§2) . In this model, observed radiation is produced by synchrotron emission of 
shock accelerated electrons. We have shown, that in the region where electrons are accelerated 
protons can be accelerated to ultra-high energy (§4.1). Acceleration to > 10 20 eV is possible 
provided that the fireball bulk Lorentz factor is large enough, T > 100, and that the magnetic 
field is close to equipartition with electrons. The former condition, T > 100, is remarkably 
similar to that inferred based on 7-ray spectra, and T ~ 300 is the "canonical" value assumed 
in the fireball model. The latter condition is commonly assumed to be valid in order to account 
for observed 7-ray emission. 

Observed UHECR flux and spectrum are consistent with a model where UHECRs are pro- 
tons accelerated to high energy in GRB fireballs, provided the efficiency with which fireball 
kinetic energy is converted to 7-rays, and therefore to accelerated electron energy, is similar to 
the efficiency with which it is converted to accelerated proton energy (§3, §4. 2). 

The GRB model for UHECR production has several unique predictions (§5). In particular, 
a critical energy is predicted to exist, 10 20 eV < E c < 3 x 10 20 eV, above which a few sources 
produce most of the UHECR flux, and the observed spectra of these sources is predicted to be 
narrow, AE/E ~ 1: the bright sources at high energy should be absent in UHECRs of much 
lower energy, since particles take longer to arrive the lower their energy. Recently, the AGASA 
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experiment reported the presence of 3 pairs of UHECRs with angular separations (within each 
pair) < 2.5°, roughly consistent with the measurement error, among a total of 36 UHECRs 
with E > 4 x 10 19 eV [jSP ]. The two highest energy AGASA events were in these pairs. Given 
the total solid angle observed by the experiment, ~ 27rsr, the probability to have found 3 
pairs by chance is ~ 3%; and, given that three pairs were found, the probability that the two 
highest energy events are among the three pairs by chance is 2.4%. Therefore, this observation 
favors the bursting source model, although more data are needed to confirm it. Testing the 
above predictions of the fireball model for UHECR production would require an exposure 10 
times larger than that of present experiments. Such increase is expected to be provided by the 
planned HiRes |3l| and Auger ^2f detectors. 



A natural consequence of proton acceleration in fireball shocks is the conversion of a large 
fraction, > 10%, of the fireball energy to a burst of ~ 10 14 eV neutrinos by photo-meson produc- 
tion (§6.1). Large area, ~ 1km 2 , high-energy neutrino telescopes, which are being constructed 
to detect cosmologically distant neutrino sources, would observe several tens of events per year 
correlated with GRBs, and test for neutrino properties (e.g. flavor oscillations, for which up- 
ward moving r's would be a unique signature, and coupling to gravity) with an accuracy many 
orders of magnitude better than is currently possible (§6.2). 
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